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Abstract: 

A manifest covariant equilibrium statistical mechanics is constructed starting with a 8N di- 
mensional extended phase space which is reduced to the 6N physical degrees of freedom using 
the Poincare-invariant constrained Hamiltonian dynamics describing the micro-dynamics of the 
system. The reduction of the extended phase space is initiated forcing the particles on energy 
shell and fixing their individual time coordinates with help of invariant time constraints. 

The Liouville equation and the equilibrium condition are formulated in respect to the scalar 
global evolution parameter which is introduced by the time fixation conditions. The applica- 
bility of the developed approach is shown for both, the perfect gas as well as the real gas. 

As a simple application the canonical partition integral of the monatomic perfect gas is 
calculated and compared with other approaches. Furthermore, thermodynamical quantities are 
derived. 

All considerations are shrinked on the classical Boltzmann gas composed of massive parti- 
cles and hence quantum effects are discarded. 
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1 Introduction and Motivation 



The motivation to investigate relativistic equilibrium statistical mechanics is given by theo- 
retical reasons as well as by practical reasons. From the theoretical point of view equilibrium 
statistical mechanics must be brought into accordance with the principles of relativity like all 
other branches of physics. From the practical point of view there might be important relativistic 
effects at (very) high temperatures. Those effects may play a major role within cosmological 
problems. Hence, equilibrium statistical mechanics is also interesting in the framework of 
general relativity and as a proper base a manifest covariant approach of a special relativistic 
equilibrium statistical mechanics is demanded. 

While the covariant generalization of kinetic equations as a major application of non- 
equilibrium statistical mechanics can be studied since many years with help of textbooks (see 
e.g. pQ) the relativistic equilibrium statistical mechanics has not reached such stable states. 
(The interest on covariant kinetic equations was mainly driven by practical reasons, i.e. by 
applications in various fields like plasma physics, astrophysics, heavy ion reactions etc.). 

Although Jiittner presented already 1911 relativistic calculations concerning the perfect 
gas [2] 2 the manifest covariant description of equilibrium statistical mechanics is still under 
consideration and a lot of attempts have been performed already to tackle the related problems. 
An overview about such problems may be studied in older papers and reviews like the ones of 
Hakim [3] and the one of Havas [3] as well as references therein. 

In this paper we mainly focus on the problem that the phase space J dT must be conceived 
in an invariant fashion demanded e.g. by the invariance of the entropy (i.e. S = klnT). Ter 
Haar and Wergeland [5] have summarized very clear the two possibilities which intuitively offer 
themselves: 

"1) Either one makes T explicitly invariant extending it from QN to 8N dimensions. That 
implies a many-time formulation of the underlying dynamics. For equilibrium theory, at any 
rate, this does not seem to be a suitable approach. 

2) Or, one retains (in the classical limit case) the usual phase space and considers the Lorentz 
transformations as canonical. This will certainly accomplish the desired invariance, but requires 
certain auxiliary re- interpretations." 

The work of Balescu et al. [U| strongly enforces possibility 2) but applying the method on 
non-equilibrium statistical mechanics only. 

The refusing statement of ter Haar and Wergeland concerning possibility 1) was certainly 
motivated by the fact that any straight forward attempt to formulate classical dynamics in a 
covariant fashion runs automatically into conflicts defined by the famous no-interaction theorem 
[7j. This theorem states that only a collection of non- interacting particles can be described 
within a Poincare invariant fashion if a straight forward generalization of the usual dynamical 
description is applied. But the conflicts with the no-interaction theorem can be avoided within 
the framework of constrained Hamiltonian dynamics. This formalism is based on an algorithm 
which was introduced by Dirac [S] and Bergmann [5j in order to express a theory based on a 
singular Lagrangian in a generalized phase space approach with help of constraints 3 . In this 
context the original 8N dimensional extended phase space given by N interacting relativistic 
particles can be reduced to 6N dimensions with the help of 2N constraints fixing the energies 
and the time coordinates of these N particles. Doing this a 6N dimensional hypersurface in the 

2 Note, Jiittners calculations are relativistic but not covariant. 

3 A general and introductional overview about this formalism can be found in monographs and textbooks 
(see HU). 
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original 8N dimensional phase space is denned on which the system is allowed to move during 
its evolution. The evolution of the system is described in respect to a global scalar evolution 
parameter which is introduced by the time constraints relating the individual time coordinates 
to this parameter. This formalism is not changing the notion of simultaneity by a change of 
the frame of reference, means one gets an invariant notion of simultaneity as well as invariant 
world lines. This method has been successfully applied already in complicated non-equilibrium 
situations namely, in order to calculate heavy ion reactions at intermediate and high energies 

mm)- 

Hence we do not agree that possibility 1) is not a suitable way and we propose in this 
paper an approach which uses the formalism of constrained Hamiltonian dynamics to describe 
the micro-dynamics of the system and which develops on this base an equilibrium statistical 
mechanics in a manifest covariant fashion. 

A similar approach has been presented already by Hsu [T3]. But this approach is based 
on the so-called 'common relativity' (developed by the same author) which is a rather special 
formulation and interpretation of relativity. In contradiction to the work of Hsu we prefer in 
this paper to apply the usual framework and interpretion of relativity. 

Another interesting approach following possibility 1) is the one of Horwitz et al. Jl] which is 
based on the many-body theory of the relativistic mechanics proposed by Stiickelberg The 
main characteristics of this method are the separation of the center of mass motion in the same 
way as in the non-relativistic theory and the allowance of the particles being off-shell resulting 
finally in a new potential corresponding to the mass degree of freedom of the relativistic system. 
However, in our opinion, the extended phase space is not reduced in a satisfying way within 
this approach: phase space integrations are only performed keeping a common time interval At 
because no specific time fixations have been incorporated in this model (see e.g. eq. (3.25) in 

El)- 

Furthermore, we would like to mention the paper of Tretyak JTB]. This author used the front 
form of dynamics to investigate relativistic equilibrium statistical mechanics. But shrinking the 
phase space ad hoc to 6iV dimensions this approach, although very elegant and usefull, can not 
be regarded to be manifest covariant. 

The presentation of our approach in this paper is organized as follows: 
In the following section 2 we give a brief overview on the constrained Hamiltonian dynamics 
in a representation which is suitable for our purpose. In section 3 we propose the formulation 
of a covariant equilibrium statistical mechanics using the framework explained in section 2. In 
section 4 we apply the developed approach to the canonical ensemble of the monatomic perfect 
gas. The obtained results are compared with other approaches. Finally, section 5 contains some 
concluding remarks. 

Some remarks on the used notation are at order: 
Throughout the paper relativistic units are used, i.e. H — c—1. The diagonal metric tensor 
has been chosen to be (1, — 1, — 1, — 1). The usual summation convention (i.e. summation of 
repeated indices in the same expression) is applied in case of Greek indices (representing tensor 
indices) which are running from to 3. Latin indices indicate the particle index. Capital 
letters are used for quantities characterizing the complete system while lower case letters are 
used for quantities related to properties of single particles. The complete set of 8iV phase space 
coordinates will be abbreviated by (g^,p M ) = (gf , q%,Pj, ...,Pj^). A contravariant 4-vector 
a M is defined as a M = (a , a). For the Poisson bracket of two dynamical phase space variables 
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and the convention 



r A B ,y- f^_^_ _ dAdB_\ 



is used. As common within constrained Hamiltonian dynamics we also distinguish between 
weak equations and strong equations, assigned by '~' and '=', respectively. While strong equa- 
tions hold throughout the extended phase space weak equations are demanded to be valid on 
the submanifold defined by the constraints. Note, weak equations should not be used before 
all interesting Poisson brackets have been worked out. 



2 Constrained Hamiltonian Dynamics 

The non-relativistic equilibrium statistical mechanics contains as an essential part the descrip- 
tion of the micro-dynamics using Hamiltonian dynamics, i.e. by classical propagation of par- 
ticles under mutual interactions " at-a-distance" . The generalization of this non-relativistic 
particle dynamics to a manifest covariant particle dynamics is not trivial because, one has to 
know how to treat action-at-a-distance in a covariant fashion. The conceptual problems in this 
field are formulated in the no interaction theorem [7j, which states that only a collection of free 
particles can fulfill the following requirements: 

• The many particle system is described in a Hamiltonian formulation with a canonical 
representation of the Poincare group. 

• The world lines are invariant (World line condition). 

• The physical coordinates are identified with the canonical coordinates. 

At least one of these requirements must be violated if one is going to allow an interaction 
between the particles which depends only on the trajectories of the particles. 

In the seventies of the last century three possibilities have been worked out to avoid this 
negative statement of the no- interaction theorem 4 : The predictive relativistic mechanics, the 
singular Lagrange formalism and the constrained Hamiltonian dynamics. All these formalisms 
have common that constraints are used to reduce the original 8iV-dimensional phase space, 
which is given by the space and time coordinates, the momenta and the energies of the N 
particles involved. The advantage of the constrained Hamiltonian dynamics is that one can 
choose these constraints more freely as in the other methods. Furthermore, one can define 2N 
constraints in a way obtaining the usual 6N dimensional phase space in the non-relativistic 
limit, containing the physical dynamical degrees of freedom. The constrained Hamiltonian 
dynamics for particles under interaction was developed by Bergmann, Anderson, Goldberg[D], 
Todorov JH], Rohrlich [T^], Komar [20], Sudarshan, Mokunda, Goldberg [21] and Samuel 
These authors picked up the idea of Dirac jH] who has realized for the first time that constraints 
are not only reducing the degrees of freedom but can also determine the dynamics. 

The model of Samuel contains as an important feature the cluster decomposition in a sat- 
isfying way for our purpose. The cluster decomposition, also called cluster separability, makes 
sure that when a system of interacting particles breaks into dynamical independent clusters, 

4 An overview about the various attempts are collected in |17j . 
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then the set of constraints has to break into these dynamical independent clusters as well. 
The disadvantage of the model of Samuel is that it is not quantizable. But this is playing 
no role for our purpose because, shrinking the considerations on the classical Boltzmann gas, 
we need only a classical Poincare invariant description of particle dynamics in the framework 
of an action-at-a-distance theory which respects the cluster separability and the principle of 
causality. 

In the following we demonstrate the concrete usage of this formalism to describe the micro- 
dynamics of the perfect gas as well as of the real gas composed of massive particles. 



2.1 Micro- Dynamics of the Perfect Gas 

Modeling the micro-dynamics of the perfect relativistic gas a system of N non-interacting 
particles can be considered. The singular Lagrangian of such systems leads directly to N 
on-shell constraints (see e.g. [23] ) 

<A = ^-(p'-O~0 ; i = l,...,N, (2) 

with being the 4-momentum of particle i and rrii being its rest mass. In contradiction to an 
earlier work and the common usage the mass-shell constraints have been modified by 
the scalar factor l/2mj. This modification provides several advantages: 

First, the transition to the non-relativistic limit is simplified and needs no further manipulations. 
Second, this factor regulates the dimensions of the mass-shell constraints to be the one of an 
energy instead of an energy squared. 

Third, and this is our major motivation, a straight forward performance of phase space integrals 
is guaranteed by this modification. 

Since these primary constraints are first class, i.e. {<fi,<Pj} ~ 0, no secondary constraints exist 
in this case. Nevertheless, the (fi have to be supplemented by N time fixation conditions which 
gauge the individual time coordinates parameterizing the world lines. Because non-interacting 
particles must not be correlated one often uses the simple gauging conditions 

Xi = q-~r^0 ; i = l,...,N, (3) 

which force the time coordinates of the particles to be equal to the global scalar evolution 
parameter r in any inertial system. These constraints are obviously not covariant. 
But our aim in this paper is a manifest covariant reduction of the extended phase space and 
hence the constraints (jHJ are not satisfying in order to treat the relativistic perfect gas. There- 
fore, we define here another set of time fixations 

X, = ^-r^0 ; i = l,...,N, (4) 

TO; 

which are covariant and also do not create any correlation among the particles under consid- 
eration. The acting of these gauging conditions is easy recognized when considering it in the 
rest frame of a particle i: In this particular frame of reference = (rrii, 0) and \i reduces to 

Xi = Qio - r ~ 0, (5) 

i.e. the time coordinates of the particles equal the global evolution parameter r in their rest 
frame. In other words, using (j3J) the proper times of the particles are synchronized to r within 
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this covariant treatment. Note, both time fixations, (J3J) and are first class, i.e. {xi, Xj} ~ 0, 
but the complete set of 2N constraints is second class since {<fi, Xj} 7^ 0. 

Because the canonical Hamiltonian derived by a Legendre transformation from the singular 
Lagrangian of a system of relativistic particles is identical vanishing j2B] the Dirac Hamiltonian 
can be constructed by a linear combination of the r-independent constraints: 



N 



H = J2^m- (6) 



This Hamiltonian generates the equations of motion in respect to r to be 



rrr X^i d( P. 




J'=l 

The Lagrange multipliers Aj can be determined with help of the consistency conditions 

f = ^ + {W«^ + EMMi}«0 ; j = l,.-.,2iV; (9) 

which guarantee the validity of the constraints during the complete evolution of the system and 
which have to be valid for the complete set of 2N constraints 

^■ = r j i0< ^ do) 

N < j < 2N. 

Introducing the invertible iV x iV-matrix C with elements = {fi,Xj} the Lagrange multi- 
pliers are uniquely determined by 

to be Aj = rrii/p® if using (j3J and Aj = nrf/p? if using (J3|. 

In the following we will use the term semi- covariant approach whenever the set of 2N 
constraints is assumed to consist of the N on shell constraints (J2J) and the simple time fixations 
because N covariant constraints are supplemented by N non-covariant constraints. Further- 
more, we will use the term full covariant approach whenever the on-shell constraints (J2J) 
are combined with the covariant time fixations (@J) to build the complete set of 2N constraints. 

Using p\ = Si + rrii the on-shell constraints reduce in the non- relativistic limit (e^ Cm;) to 

<Pi = £ i ~ 7T~ « 0- (12) 
2mi 

Furthermore, the covariant time fixations (p£J) reduce to the simple time fixations © in the 
non-relativistic limit and the global evolution parameter r is identified with the absolute time 
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t. Consequently, the usual non-relativistic Hamiltonian dynamics is obtained as generated by 
the non-relativistic Hamiltonian 

N 2 

*0-> = ££-. (13) 

i=l 



2.2 Micro-Dynamics of the Real Gas 

In contradiction to ideal systems like the perfect gas real systems can only be understood if the 
intermolecular interactions are taken into account when describing the related micro-dynamics. 
Applying the formalism of the constrained Hamiltonian dynamics to the micro-dynamics of 
the real gas we therefore have to consider a system of N particles in mutual interaction. In 
this case we prefer to use a set of constraints quite similar to the constraints which have 
been successfully applied already within non-equilibrium situations when simulating heavy ion 
reactions (see JI],P2])- The first N constraints are chosen as on-shell constraints 

Vi=2^:(p?-m?-y i )«0 ; i = l,...,N. (14) 

Note, as already discussed in section 2.1, the on-shell constraints have been modified by the 
scalar factor l/2mj. Regarding the potential part Vi this choice of the constraints requires that 
V{ should be a Lorentz scalar and therefore a function of Lorentz scalars. Defining a system 
with mutual two-body interactions, V should be given by a sum of these two-body interactions 

Vi = Y.ViMhi)- (15) 

Using equation (JT5J) we assume that the two-body interactions depend only on the Lorentz 
invariant squared transverse distance 



2 2 (lijPijn) 

Zra = % -2— > ( 16 ) 



with = q^ — qj being the four dimensional distance and = p^+Pj the sum of the momenta 
of the two interacting particles i and j. This particular choice has two advantages: 
First, in the common CMS of two interacting particles CMSy q\^ reduces to — q?- which is 
quite natural since only in the CMSjj particle % and particle j can be treated on the same 
footing. In addition, in the non-relativistic limit the potential part Vi can be connected with 
the corresponding potential V, used in the non-relativistic theory: Using p° = + rrii and 
regarding the non-relativistic limit (e, uij) the on-shell constraints reduce to 

d 2 V 

Vi = - ^ - ip~ « 0. (17) 
ivrii zrrii 

The comparison of (JT7j) with the usual expression of the non-relativistic energy suggests 

Vi = 2miVi, (18) 
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and the two-body potential acting among two particles i and j felt by particle i reads in the 
common CMSy 

Vfc(-qJ) = 2m i V ij (-c£ j ), (19) 

i.e. a non-relativistic two-body potential can be generalized to be used in the relativistic theory. 
For instance, the widely used Lennard- Jones potential can be generalized replacing the squared 
distance qf ■ by the Lorentz invariant squared transverse distance gf^ ■ obtaining 



v£ J \ql ij ) = 2m i K 




with the two 'tuning parameters' k and a having the dimensions of an energy and a length, 
respectively. Note, due to the second term of the right hand side of equation (jTBj) this generalized 
interaction is slightly implicit momentum dependent. Because of this term, which gives the 
longitudinal squared distance, the interaction depends not only on the distance of the two 
interacting particles but also on the direction of their center of mass motion in the rest frame 
of the real gas. 

Since the on-shell constraints alone do not specify the world lines one needs additional N 
constraints which are fixing the relative times of the particles. Like in the model of Samuel [22] 
this second set of N constraints should fulfill the following conditions: 

1) N—l of them should be Poincare invariant in order to fulfill the required world line invariance. 

2) Causality has to be respected. 

3) One has to be consistent with cluster separability. 

4) A global evolution parameter must be defined. 

These conditions can be fulfilled by defining the following set of N time fixations 

Xi = — ; z = l,...,(JV-l) (21) 

771/4 



Xn = P^Q^ - r « (22) 

with = P^/Vp^, = J2iPii Q 11 = an d the dimensionless scalar weighting 

function 

c^ = — ^efr^ (23) 

with a being a characteristic interaction distance of the modelled interaction. Note, the scalar 
factor 1/rm regulates the dimension of the time fixations to be the one of a length. The 
conditions ()21|) are motivated by studies in the framework of Singular Lagrangians [21]. Using 
this methods one gets up to the weighting functions Uij the same conditions as secondary 
constraints and the on-shell conditions as primary constraints in a natural way. The important 
fact for using the expression (J23)) as weighting function is that this scalar function respects the 
principle of causality, while the ones used in the Singular Lagrangian theories and in the model 
of Samuel can violate this important physical restriction (see |llj). 

The constraints (|21|) take care that the times of interacting particles are not dispersed to 
much in their common CMS^. Furthermore, they specify the dynamics by fixing the times at 
which the forces have to be calculated but they do not specify the global evolution parameter r. 
This evolution parameter must defined to be determined dynamically since the no-interaction 
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theorem can only be avoided in this way and is defined by the gauging condition (J22J) in a way, 
that the individual times are increasing with increasing r. 

In the CMS of all particles involved r is simply given by the average of all time coordinates 

r=P*Q, 25 lf>°. (24) 

1=1 

This set of 2N constraints given by (|14 p .(|2i p and (|22j) determines covariant world lines 
parameterized by the initial data at equal r. 

A linear combination of the 2N — 1 r-independent constraints 

'' 0< *- N (25) 
(Xi-N ;N<i<2N-l 

constructs the Dirac Hamiltonian 

2N-1 



H=Y, X ^ ( 26 ) 



i=l 



which generates the equations of motion in respect to r to be 

cfgf 2 v~^ 1 dip 



i=i 



The Lagrange multipliers Aj can be determined with help of the consistency conditions 



2JV-1 



+ + X) Ai{^,^}~0 ; i = l,..,2iV; (29) 



dr dr dr 

i=i 

and are finally given by 

A, » "^C- (30) 

with C being the invertible 2N x 2iV-matrix with elements = 

The whole formalism has a well defined non-relativistic limit as has been shown by Sorge 
[TT] . In leading order of 1/c one gets for the first 2N — 1 constraints 

W = ei-7T--7r-™Q ; i = l,-..,iV (31) 
2mj 2mj 

Xi = qi-q° N ~Q ; i = i,...,JV-i (32) 

which generate the same dynamics as the non-relativistic many body Hamiltonian 

N / 2 \ 

(33) 



/ 2 \ 

^" = E fe + y) • 
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3 Covariant Statistical Mechanics in the Framework of 
Constrained Hamiltonian Dynamics 

Regarding the constrained Hamiltonian dynamics as a proper formalism to describe a system 
of interacting particles in a covariant fashion we are to develop a manifest covariant statistical 
mechanics on its base. We start this development considering the extended phase space. In 
case of a system involving N particles this invariant phase space has 8N dimensions given by 
the N 4-momenta pf and the N 4-position- vectors gf (i = 1, .., N). 

The phase space distribution function D(q^ , p M , r) representing the statistical ensemble is as- 
sumed to depend on the complete set of 8N phase space coordinates and, in general, on the 
global evolution parameter r which is introduced by certain time fixation conditions (compare 
section 2). Consequently, the Liouville equation and the equilibrium conditions have to be 
discussed in respect to r. 

The Liouville equation is thus given by 

dr or 
The micro-dynamics of the system is generated by the Hamiltonian 

M 

H = J2 X ^ (35) 
i=i 

where M = N in case of the perfect gas and M = 2N — 1 in case of the real gas using the 
constraints as introduced in section 2. 
Hence we get 

ao y ,T) »EMflK.PV),ft}. (36) 

i=l 

The reduction of the invariant phase space volume element of the extended phase space 

TV 

dr e = Y[d 4 qi d%. (37) 

i=l 

is performed incorporating all 2N constraints by means of 5-functions, i.e. 

N 

dT = H d%d% 9(p°)5(^)5( Xi ). (38) 

i=l 

The Heaviside-G-function has been added in order to shrink on particles with positive energy 
which is usefull if classical particles are considered only. 

Equilibrium distribution functions are r independent solutions of ()34|) and hence have to 
fulfill the condition 

AI 

^Ai{D(^,^),Vi}«0. (39) 

t=i 

Using (|3*Hj) the ensemble average (i.e. the classical expectation value) of a measurable quantity 
A(p tM , g M ) at equilibrium is given by 

<A>= J UL d%d% Q(p°)5(<PiMXi)A(p», q»)D{q^ p") 
IIllid 4 qid%Q(p^5(<p i )S(Xi)D(q^p^ 
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We mention two simple (and widely used) possibilities to construct solutions of the Liouville 
equation fulfilling the equilibrium condition: 



• Constant phase space distributions D = const. 

• Phase space distributions depending on quantities Xj which are conserved during the 
evolution in respect to r, i.e. dXj/dr = 0. 

Assuming Xj ^ Xj{r) the condition 



has to be valid in case of such conserved quantities. 

Beside the trivial uniform ensemble the microcanonical ensemble can be regarded to be the 
most typical example for the former case while the canonical as well as the grand canonical 
ensembles are examples of the latter case. 

In the non-relativistic theory the energy is playing a prominent role as the most important 
conserved quantity. In a relativistic generalization the energy has to be replaced by the 4- 
momentum P M or even by the energy-momentum tensor T^ v of the system. A fruitfull alternate 
is the invariant energy E = U^P^ 1 with being the 4-velocity of the system. In the comoving 
frame of reference (i.e the restframe of the system) we get E = UqP° = P°, i.e. the invariant 
energy is nothing but the total (internal) energy of the system measured in its restframe. 
We will use this invariant energy when generalizing the concepts of non-relativistic statistical 
mechanics. 

In the following we consider briefly the main equilibrium distributions: 
i) Microcanonical Ensemble 

Considering the classical (i.e. non quantum mechanical) case only and applying our concept of 
the invariant energy E we can write the distribution function of the microcanonical ensemble 
using the ^-function (compare e.g. [25] ) as 



with Eq being the invariant energy of the closed system. The 'microcanonical partition integral 
Zm (which is nothing but the phase space volume T) has to be determined via 



M 




(41) 



i=l 




(42) 




(43) 



The factor 1/N\ respects the fact that the particles can not be distinguished. 



ii) Canonical Ensemble 



Considering the canonical ensemble the condition 



M 



M 




(44) 



i=i 



i=i 
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is to be proved in order to use the ansatz 

D c (E) = ±e-^ (45) 
for the canonical distribution function with the canonical partition integral 

N 

Zc = m J n e (p?)%*)%*) e_/3c/MPM - ( 4e ) 



iii) Grand Canonical Ensemble 

The grand canonical ensemble is characterized by another conserved quantity beside the energy 
namely, the particle number N. Hence, generalizing the canonical ensemble accordingly leads 
to 

D G (E,N) = ^e-^-« N (47) 
for the grand canonical distribution function with the grand canonical partition integral 

N 

Z c = Jn U e(p?)tf( W )*( X i) e-^^-^. (48) 

' ^ i=i 



iv) Related Thermodynamics 

Let us briefly discuss the corresponding thermodynamics generated by our approach. Like in 
the next section, we shrink on the canonical ensemble. The formal properties of this thermody- 
namics are dominated by the introduction of the invariant energy E. As in the non-relativistic 
theory the free energy can be determined with help of the canonical partition integral Zc by 

F = -kT\nZ c . (49) 

Like the invariant energy E the free energy F (as all other thermodynamical potentials) is 
defined to be a scalar. Furthermore, the temperature T introduced by the parameter (3 = 1/ kT 
(k being the Boltzmann constant) is a scalar quantity and the introduction of a heat vector 
[2*o] is superfluous 5 . Considering the usual relations defining other macroscopic state variables 
like the entropy 

and the pressure 

'"(§& 

5 Of course, formally we get the same exponential factor (see eq. I@SJ) but the 4-velocity of the system U ^ 
is now defining the invariant energy E = U^P^ and not used to define the heat vector M = ^ like in other 
approaches |20| . One can argue about this interpretation but the advantage defining an invariant energy is that 
this scalar quantity of the exponent can be demanded to be conserved during evolution in respect to r, i.e. to 
fulfill equation 
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we stress that both, the macroscopic state variables as well as the thermodynamical potentials 
are defined to be scalar quantities in our approach. Hence, an artifical extension of those 
quantities to higher rank tensors as often done when formulating relativistic thermodynamics 
(compare e.g. is not needed. 
The average energy is given by 

<-> - m 

and the specific heat can be determined with help of it: 



dT 



v 



Looking to the grand canonical ensemble we note that the chemical potential M introduced 
by the parameter a = (3M = M/kT is a scalar as well in the approach considered. 



4 The Canonical Ensemble of the Monatomic Perfect 
Gas 

In this section we want to apply the formalism developed in the previous sections. Unfortu- 
nately, concrete calculations (e.g. determining the canonical partition integral via equation 
(jjBj) ) are note easy to be performed due to the structure of the phase space volume element 6 . 
Hence, as an example, we just apply the formalism to the canonical ensemble of the monatomic 
perfect gas. In order to determine the canonical partition integral using formula 1)46)1 we have 
to prove the validity of condition (|4"4"j) : 

By definition, the perfect gas is a system of non-interacting particles and as a consequence the 
Hamiltonian contains the simple mass shell constraints (JSJ). Hence, the partial derivatives of ifi 
in respect to g JM are vanishing and condition (|44|) is fulfilled in case of the perfect gas. 

Unfortunately, no analytical solution of the canonical partition integral ()46j) can be achieved 
within the full covariant approach. Hence, numerical solutions have been worked out as will be 
explained in section 4.2. 

But before discussing these numerical results let us first present another approach, namely the 
semi-covariant approach which provides analytical results and hence can be easy compared with 
other approaches known from the literature, e.g. the non-covariant Jiittner approach 



4.1 Semi-covariant Approach 

In the semi-covariant approach the simple time fixation conditions are used and the partition 
integral reads 

z ° = m / n d ^ (^-^ 2 - ™ 2 )) w - r ) e ~" u ^- ^ 

6 Especially, as in the non-relativistic theory, calculations considering the real gas turn out to be complicated 
depending on the underlying model of intermolecular interactions. 
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Using the properties 



and 



of the 5-f unction we get 



5(ax) = y-r6(x) 



5(x 2 — a 2 ) = — [5(x — a) — 5(x + a)] 



Zc= w\ /n ^ (rf - H + m 5 w - r ) e ~^ pM - 



(55) 
(56) 

(57) 



For the further concrete calculation we assume a comoving frame of reference (i.e. the restframe 
of the gas). In this frame of reference we get 



JV 



N 



(58) 



i=l i=l 

and integrating the O-components gives according the 5-functions 

TV 



am 



ft 



JV 



mi 



8=1 



2 



(59) 



where V represents the volume measured in the restframe of the gas. The integral factorizes 
using spherical coordinates in momentum space and performing the angle integration gives 



A? 



V 
AH 



47rm / dp 



N 



(60) 



with p = yp^. Note, the particle index becomes superfluous considering a monatomic gas 
of particles with mass mj = m. The integral in this expression can be solved applying the 
transformation p = msinh£ (see e.g. |27j ) 



dp 



m 

[3m 



(61) 



with K\ being the modified Bessel function of first order. 

Finally, the canonical partition integral within this semi-covariant approach is given by 



V 
AH 



(3m 



N 



(62) 



14 



Now we are able to list the analytical results of 3 different approaches: 



; semi-covariant 

; Jiittner approach (53) 
; non-relativistic (see e.g. L 28j). 

Comparing (J52*|l with the result of Jiittner we realize that the only difference is the order of 
the modified Bessel function occurring in Zc'- Jiittner 's approach contains K 2 {l3m) instead of 
K x {(3m). 

We note further that this result is naturally quite similar to the one derived by Horwitz et al. 
(compare eq. (3.25) in [14J, especially regarding the functional dependence on (3. But we stress 
that the phase space reduction regarding the time constraints in our semi-covariant approach 
is already more satisfying and as a consequence, in contradiction to [Hj, no commen time 
interval is remaining in the final result when applying the formalism of constrained Hamiltonian 
dynamics. 



yN 


Anm 3 , 


m 


(3m 1 


yN 




m 


/3m 


yN 


1 2-ivm 


m 





1 N 



N 



3/2 



1 N 



4.2 Full covariant Approach 

Evaluating the partition integral within the full covariant approach we apply the covariant time 
constraints (J3J as specified in section 2.1 for a system of non-interacting particles and equation 
reads 

Z ° = M / A d ^ - 5 - r) e-^. (64) 

Using the properties (J55J) and (jSUjl of the 5-function we get 

z ° = m / n (p" - vp? + m i) 6 (9? - - mr )) e ~" u ^- ( 6s ) 

Until yet we have kept strictly the covariant fashion during the calculation. For the fur- 
ther concrete calculation the fixation of a defined frame of reference is demanded. Hence, for 
simplicity, we choose the restframe of the gas to determine the invariant energy E as already 
applied within the semi-covariant approach (see eq. ([58))) and get integrating the O-components 
according the ^-functions 




The integral factorizes using spherical coordinates in momentum space and performing the 
angle integration gives 

yN 

z c = 

c m 



Airm 



dp 



p2 _|_ m 2 



(67) 
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with p = yp^ Note, as in equation (|60jh m; = m has been used again. 

Unfortunately, no analytical solution of the integral in equation ()67|) is available. Hence, 
numerical solutions have been worked out for various monatomic gases. 

In order to compare this numerical results obtained by the full covariant approach with the 
results of other approaches as listed in equation (|63|) we write the canonical partition integral 
of the monatomic gas in a generic representation: 



c = — (e~^Y) N . (68) 



yN 

iV! 

The 'relativistic renormalization factor' e _/3m has been extracted explicitly in (J68|) for two 
reasons: 

1) This factor does not occur in the non-relativistic partition integral. 

2) The numerical integration is simplified if this factor is extracted explicitly. 

The non-analytical integral occurring in the full covariant approach (see eq. (JfiTfO has been 
determined numerically using the Gauss-Laguerre method of 15 th order. The results of the two 
other relativistic approaches (i.e. the Jiittner approach and the semi-covariant approach) have 
been determined using the analytical expressions (i.e. evaluating the related modified Bessel 
functions K\ and K-i) as well as using the same numerical integration method in order to justify 
the accuracy of the applied numerical method. 

Figure 1 summarizes the results of calculations for different monatomic perfect gases name- 
ley, hydrogen, helium, neon and argon. The comparison is focused on the quantity Y which 
determines the canonical partition integral via equation (pjH|) . As expected, the qualitative 
behavior is the same for all gases. 

At low temperatures (non-relativistic limit) all approaches give the same result as it should 
be. Even at higher temperatures (T < 10 12 K) 7 no difference is visible among the various 
calculations. At extremely high temperatures (T > 10 12 i^) relativistic effects are dominating 
the behavior. Compared to the non-relativistic result the non-covariant Jiittner approach is 
showing a strong enhancement of the quantity Y in this regime. The semi-covariant approach 
reduces this enhancement and the full covariant approach turns the enhancement to a reduction 
even. Therefore, we conclude that a controlled manifest covariant reduction of the extended 
phase space reduces the value of the quantity Y and hence of the canonical partition integral. 

In order to give an overview on the related thermodynamics of the monatomic perfect gas 
we finally list in dependence of the quantity Y the free energy 



the entropy 



F = -kT\uZ< • = -M,T ( In— + \nY - (3m + 1 ) . (09) 



\dT J v dT 



( V d\\iY \ 

= M^_ + ] n y-0_ + lj (70) 

the pressure 

fdF\ , d\nZ c NkT , . 

P — — I — — I = kT—-^- = , 71 

\dV T dV V y ! 



7 E.g. T = 1Q 12 K corresponds roughly to /3m ~ 10 in case of hydrogen. 
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Figure 1: Quantity Y as a function of (3m for various approaches and different monatomic 
perfect gases. 

the average energy 
and the specific heat 

cv= ( a -m =^^. (73) 



dT J v T df3 2 

Note, since the pressure turns out to be independent of Y all approaches give the same 
equation of state, namely the one of the non-relativistic perfect gas: 

PV = NkT. (74) 

Since the main differences among the various approaches are showing up at very high tem- 
peratures it is very usefull to consider explicitly the 



4.3 Ultra-relativistic Limit 

Since in the ultra-relativistic limit the rest mass is neglected in respect to the momentum and in 
respect to the energy it can be evaluated using yp? + mf — > |pj| when deriving the integrals 
in eqs. (16 7|) and (JoTij) . Hence, it is even possible to derive analytical results in case of the full 
covariant approach since the partition integral simplifies drastically and we get in this limit 
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(p = |p|) for the various approaches: 



[4vrm 2 f^dpe-Pp] 
[47im f™dppe~ /3p ~\ 
[4n f-dpp^e-^Y 



N 



full covariant 



y(ur) 



N 



semi-covariant 



(75) 



N 



non-covariant (Jiittner) 



The integrals in (|7H)l are easily determined with help of the T-function according 



/ d(Ce 



-aC - 



r(n+l) 



(76) 



Table 1 summarizes these results as well as some derived thermodynamical quantities in case 
of the ultra-relativistic limit. Progressing from the non-covariant approach towards the full 
covariant approach via the semi-covariant approach at a fixed (3 the canonical partition inte- 
gral is decreasing as already mentioned when discussing figure 1. This behavior as well as the 
derivatives in respect to (3 are fixing the differences in the thermodynamical quantities, e.g. the 
entropy is reducing when progressing towards the full covariant approach. Especially, we want 
to focus on the average energy: 

While Jiittners non-covariant calculations at the ultra-relativistic limit lead to 3NkT the semi- 
covariant approach (as well as the method used in [Tlj) leads to 2NkT and the full covariant 
treatment gives NkT. Thus, a controlled manifest covariant reduction of the extended phase 
space by N on-shell constraints reduces the average energy by a factor NkT and the final re- 
duction by iV time fixations by another factor NkT. Consequently, the specific heat is showing 
similar differences. As already mentioned by Horowitz et al. JH] no empirical evidence distin- 
guishing between these results is available at the present time. 
The following should be noted: 

Usually, the non-covariant Jiittner approach in the ultra-relativistic limit is regarded to rep- 
resent the classical photon gas. This approximate interpretation is not possible within our 
approach because the applied classical constrained Hamiltonian dynamics of massive particles 
can not be used in the ultra-relativistic limit to approximate massless classical particle dynam- 
ics. This fact is showing up in the results and can be recognized inspecting Zq (see table 1): 

In contradiction to the Jiittner approach Zq remains mass dependant in the ultra-relativistic 
limit in case of the semi-covariant and the full covariant approach. 



4.4 Non-relativistic Limit 

Finally, we examine the non-relativistic limit of our approach. Since the covariant time fixations 
P| reduce to the simple time fixations (J3J) in this limit the canonical partition integral within 
both approaches, semi-covariant and full covariant, coincide according eq. (fT2J) to be 



Special care has to be taken regarding the invariant energy E. In order to get a fair comparison 
in the non-relativistic limit the rest masses of the particles have to be subtracted from their 




(77) 
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thermodynamical 
quantity 


full covariant 


semi-covariant 


non-covariant 
(Jiittner approach) 


y{ur) 


r i N 

V N 4vrm 3 

AT! [ (/3m) 


V N \ 4irm 3 ] N 
N\ [(/3m) 2 J 


r iiV 

V N 8wm 3 


F = —kTlnZc 


a t i rn n V i i / a 3 \ 

— N kT[mjj + m(47rm cS ) 
-ln(/3m) + 1] 


A r 7 rn n 1/ i 1 /i 3 \ 

— A/ fcT[m-^ + m(47r?7r) 
-21n(/3m) + 1] 


A77 rn l~1 \f i 1 /o '"i \ 

— NkT[mjf + ln(87rm d ) 
-31n(/3m) + 1] 


S = klnZc 
-±dlnZ c /d(3 


iVfcpn-^ + ln(47rm 3 ) 
-ln(/3m) + 2] 


Nk\\n% + ln(47rm 3 ) 
-21n(/3m) + 3] 


A^A;[ln^ + ln(87rm 3 ) 
-31n(/3m) + 4] 


(E) = -d\nZ c /d(3 


NkT 


2iVA;T 


3NA;T 


c v = d(E)/dT 


Nk 


2NA; 


3NA; 


p = -dF/dV 


NkT/V 







Table 1: Thermodynamical quantities for the perfect gas in the ultra- relativistic limit 
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energies getting in the comoving frame of reference 



N 



i 

Using further p° = E{ + m 8 equation (J77)l simplifies to 



(7* 



7 (nr) 



^ / 



JV 



1=1 



~/vT 



2nm z 

(3m 



2 \ 3/2' 



iV 



(79) 



i.e. the non-relativistic canonical partition integral. 

Note, the same result can be achieved starting with the analytical result as obtained by the 
semi-covariant approach. Supplementing eq. (|62j) by a factor e^ Nm which covers the effect of 
the reduction of the energies by the rest masses (compare eq. (fTKJ) ) equation (JH2|) transforms 
to 



N 



V 



(3Nm 



(3m 



KtiPm) 



N 



(80) 



In the non-relativistic limit (m ^> kT) the asymptotic form of the modified Bessel function 
((3m — > oo) 

An 2 



K n ((3m) 



71 



-f3m 



2 (3m 



1 + 



l(3m 



+ 



(81) 



can be applied. Since K\ and K 2 coincide in the leading order the semi-covariant approach 
and the non-covariant Jiittner approach provide the same result in the non-relativistic limit 
namely, the non-relativistic canonical partition integral (J75j) within this treatment. 
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5 Summary and Conclusion 



A covariant equilibrium statistical mechanics has been formulated on the base of constrained 
Hamiltonian dynamics. The usage of this formalism to describe the micro-dynamics of the 
ensembles guarantees the manifest covariance of the developed approach. 

Unfortunately, concrete calculations are not easy performed within this approach. As a 
simple application the canonical partition integral of the monatomic perfect gas has been eval- 
uated numerically. Relativistic effects at very high temperatures have been observed: 
The comparison of the obtained results with the results of the non-relativistic theory show a 
decreasing partition integral at very high temperatures. In contradiction to this findings other 
relativistic approaches like the non-covariant calculations of Jiittner are showing an increased 
partition integral at these temperatures. It has been demonstrated that the full covariant 
treatment of the phase space is responsible for this significant difference. Considering the ultra- 
relativistic limit the same behavior is visible on a full analytical base. In the non-relativistic 
limit the results coincide with the results obtained within the non-relativistic theory. 

As in the non-relativistic theory the link to thermodynamics has been performed via the free 
energy. All thermodynamical potentials (like the internal energy, the free energy, ...) as well as 
the thermodynamical state variables (like temperature, entropy, pressure, ...) have been defined 
to be Lorentz scalars. This was achieved introducing the invariant internal energy E = U^P^ 
in the canonical distribution function. Regarding the ultra-relativistic limit the entropy, the 
average energy and the specific heat are decreased in comparison to other approaches in this 
limit . 

Unfortunately, the temperature regime which shows significant differences among the vari- 
ous approaches (T > 10 12 i^) can not be reached in experiments but may exist in cosmological 
events. Hence, it may be interesting to use the presented special relativistic approach as a base 
to develop a general relativistic equilibrium statistical mechanics being manifest covariant from 
first principle. 
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